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ABSTRACT 

Analysis of the Pangu N-body simulation validates that bulk flow of halos follows Maxwellian distribution 
of which variance is consistent with prediction of linear perturbation theory of structure formation. We propose 
that consistency between observed bulk velocity and theories shall be examined at the effective scale as radius 
of spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the 
sample window does. Then we compared some recently estimated bulk flows from observational samples 
with prediction of the ACDM model we used, some results deviate the expectation at level of ~ 3a but the 
tension is not as severe as previously claimed. We disclose that bulk flow is weakly correlated with dipole of 
internal mass distribution, alignment angle between mass dipole and bulk flow has broad distribution but is 
peaked at ~ 30-50°, meanwhile bulk flow shows little dependence on mass of halos used for estimation. In 
the simulation of box size l/r'Gpc, for a cell of radius 100 _1 Mpc the maximal bulk velocity is > 500kms _1 , 
dipoles of environmental mass outside the cell are not tightly aligned with the bulk flow, instead are located 
randomly around it with separation angles ~ 20-40°. In the cell showing largest bulk velocity there are slightly 
smaller number of low mass halos, however halos inside are clustered more strongly at scales > 20/i _l Mpc, 
which might be a significant feature since the correlation between bulk flow and halo clustering actually grows 
into notable beyond such scales. 

Subject headings: galaxies:halos — large-scale structure of universe — methods: statistical 



1. INTRODUCTION 

Bulk flow refers to the apparent coherent peculiar motion 
of galaxies and galaxy clusters in a considerably large vol- 
ume around us. In practice there are several ways to estimate 
bulk flows from various observation resources, such as galaxy 



catalo gues from peculiar velocity survey (e.g. iFeldman et al 
2010), compiled Type la supernovae data (e.g. iDai et al 



201 II) . galaxy clusters in combinatio n with cosmic microwav e 
background (CMB) observation (e.g. Kashli nsky et al.ll2010l) . 
In recent there emerge also interesting new methods base d on 
galaxy two- point correlation function dSong et al.ll20TTb and 
galaxy light dNusser et ai]|201 HI Abate & FeldmarfeoA 

A nalysis of the spiral g alaxy catalogue of the SFI++ sur- 
vey (ISpringob et al.l l2007) shows that within a top-hat spher- 
ical window of radius 40/r'Mpc bulk flow is 338 ± 38kms _I 
towards Galactic (l,b) = (276°, 14°) with 3° la uncertainty, 
and then 257 ±44kms _I towards (l,b) = (27 9°, 10°) with 6° 
error within window of radius 100/r'Mpc (Nus ser & Davis! 
2011b. The measure ment is in agreement with analysis of 



Sanda geet al.l (|2010) to data consisting of s upernovae, se- 
lected nearby galaxies and galaxy clusters. Feldm arTet al.l 
(2010) construct a composite catalogue of galaxies with pe- 
culiar velocities measured in different surveys, including the 
SFI++, estimate that bulk flow within a Gaussian window of 
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50/r'Mpc is 416 ± 7810ns" 1 in directio n of (l,b) = (282 ± 
1 1 ° , 6 ± 6°) (also see IWatkins et al]|2009h . 

Employing peculiar velocities of supernovae is another vi- 
able route to detect bulk flow, though such sam ples usu- 
ally ar e very sparse and prone to Malmquist bias. IDai et all 
(1201 ll) fitted a bul k flow to the Union2 supernovae catalogue 
(lAmanullah et al J 120101) of 188!j^kms _1 in the direction of 

(l,b) = (290^f ,20!|°) within redshift z < 0.05, while no 
significant bulk flow is det e cted f rom data at z > 0.05; with 
the same data IColin et all d201 ll) obtain similar estimation 
but with a higher median amplitude 250^JgQkms _1 . How- 
ever, using a differe nt supernovae data se t within redshift shell 
z = (0.0043 , 0.028). IWeyant et al.l (1201 ll) estimates that the lo- 
cal flow is 538 ± 86kms _1 pointing to (/, b) = (258 ± 10°, 36 ± 
ll°),or446±101kms- 1 towards (l,b) = (273 ± 11°,46±8°) 
if a different technique employed, which is in agreement with 
the dipole of CMB ( l b) = (26 3 99 ± 0.14° 48.26 ± 0.03°) 
(Jarosi k et all 120111) . Uhaet al.l (120071) and Haugb0ll eetaLl 
(1200717 have the same figures with analogous uncertainty. 

Availability of modern galaxy peculiar velocity data is 
limited to our local Universe, bulk flow at higher redshift, 
sometime dubbed dark flow, is mainly explored through 
kinetic Sunyaev-Zerdovich ( kSZ) effect of ga l axy c lusters 
dSunvaev & Zeldovich|[T980T) . iKashlinskv eTall (1201 ll) com- 
pute kSZ signals of 771 x-ray clusters in the 7 yr WMAP 
CMB map, conclude that the flow at z < 0.16 is directing to 

(/,6) = (278±18°,2.5 ±15°) and then (/,£) = (283 ±19°,20± 
15°) if z < 0.25. They further argue that the flow at these 

depths shall reach magni tude of ~ 1 000km s" 1 a ccording to 
an earlier investigation in Kashlins ky et al.l (12 010). But a con- 
flicting assertion is derived by Osb orne et al.l d201 ll) from the 
same CMB map in conjunction with 736 ROSAT observed 
clusters , which tells that there is no significant detection of 
kSZ effects at low multipoles and basically deny the existence 
of bulk flow, meanwhile in some cases the thermal Sunyaev- 
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Zel'dovich effect might induce a dipole that could be easily 
misunderstood as bulk flow of ~ 2000-4000kms _1 . 

Bulk flow is a topi c of long term interest s to observa- 
tional cosmology (see IStrauss & Willicklll995l for a review 
of early w orks), to which ther e are special surveys dedi- 
cated (e.g. iCourtois et ail 1201 11) . However, as we can see, 
no consensus on amplitude, direction, convergence depth of 
bulk flows has been achieved currently to reconcile differ- 
ent measurements. Nonetheless, some authors have argued 
that the amplitude of their measured bulk flow is too strong 
over such large scales, presenting a challenge to the standard 
ACDM model or at l east the one of 5yr WMAP parameters 
(Watk ins et all [20091: iKashlinskv et alj l2010t iFeldman et all 
l2010t iMacaula v et al.112.01 lh . Not strangely there have been 
cosmological models of different flavors constructed to ex- 
plain such anomaly (e.g. iMersini-Houghton & Holmar3l2009l; 
lAfshordi et aT]|200^IWyman & Khourvll2010l) . 

Expectation about bulk flow in ACDM universe is gener- 
ally calculated with linear perturbation theory of large scale 
structure in that at large scales its accuracy is believed suf- 
ficient. Often the 1-D velocity variance of dark matter are 
quoted to compare with measurements, however we need to 
address here that it is the rms velocity that should be used 
instead of the 1-D rms velocity. Linear theory predicts that 
the rms bulk velocity is typically ~ 300kms _1 with la cos- 
mic variance approximately 200km s -1 for a top-hat window 
of radius 60/i _I Mpc (Mak et aLll201 11) . A concern is that non- 
linearity might no t be negligible even at very large scales (e.g. 
Scoccimarro 2004), which could act as systematical bias to 
the conclusion about the consistency between model and data. 

Although linear theory can predict the possibility of ob- 
serving a bulk flow of particular amplitude at certain scale, 
several key problems yet can not be easily tackled analyti- 
cally, e.g. internal properties of the volume demonstrating 
large bulk blow. Practically halo catalogues from N-body 
simulation in large box with sufficient mass resolution are 
best suited for such task. The reason of focusing on halos 
instead of dark matter is that observational objects used to de- 
termine bulk flow are galaxies and galaxy clusters which are 
residing in halos, and in practice the strongly non-linearity of 
peculiar velocities of galaxies is largely filtered out so that 
what contribute to bulk flow es timation is mainly th e motion 
following the i r host halo s (e.g. IWatkins et al.l 120091) . In fact 
iBahcall et ail dl994 and Moscard ini et al.l (Il996l) "have per- 
formed analysis of mock halo catalogues and obtained useful 
results, but their simulations are either of very low mass res- 
olution or based on compromised simulation method. In this 
paper we will demonstrate our analysis of the velocity field 
of halos resolved from a dark matter only ACDM simulation 
in l/r'Gpc box with 3072 3 particles. The large volume and 
high mass resolution of our simulation enables investigating 
halo behaviors in detail over broad dynamic ranges supersed- 
ing previous works. 

In section 2 definition and basic theoretical prediction of 
bulk flow is introduced, which is followed by section 3 pre- 
senting measurements of bulk flow of randomly placed cells 
in simulation. Section 4 is devoted to analysis of special re- 
gions showing extraordinarily large bulk flow velocity. Sum- 
mary and discussion is in the last section. 

2. BULK FLOW: DEFINITION AND LINEAR THEORY 

Placing a window of characteristic scale R randomly in the 
sample space, if N objects (galaxies or halos, in this work just 



the latter) are enclosed, bulk flow of the particular volume 
indicated by the particular object is 



Wi 
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in which v, is the peculiar velocity of the ith object and w,- is 
the weight assigned. Practically the weights could be origi- 
nated from radial selection function, angular selection func- 
tion, luminosity, mass and etc. In simulation we do not have 
those observational systematical biases, we will discuss just 
two weighting schemes: window function and halo mass. 

Weighting by window function only is pure geometrical, 
which can be easily modeled. In continuous limit, 



\ = y(g)W(R) = 



l 



(2tt) 3 



\(k)W(kR)d 3 k , 



(2) 



where W(kR) is the Fourier transform of the window func- 
tion W(R). The mostly common are top-hat window 
Wth = 3(sin kR-kR cos kR)/(kR) 3 and Gaussian window Wc = 
exp(-k 2 R 2 /2). There are not much differences applying top- 
hat or Gaussian window in principle for bulk flow measure- 
ment if effective volume and depth has been taken care of. 
Sometimes bulk flow estimation is provided by objects within 
a spherical shell defined by two radius R\ > Rq, it is easy 
to see that Eq. [2] applies with window function Ws = (R\Wi - 

RlW )/(R]-Rl). 

Probability distribution function (PDF) of V could be ex- 
pressed as 

p(Y)d\ = p(V)dVp(n v )dn v (3) 

where hy denotes the unit vector in direction of the bulk veloc- 
ity. Isotropic assumption leads to p{hy) = 1 /A-tt and (ny) = 0, 
which ensures that (V) = J Vp(V)dV J n v p(n v )dn v = (V) ■ 
(hy) = 0. Integration over the angular part of Eq. [3] yields 
the PDF of the amplitude of bulk flow 



p(V)dV = 



p(Y)dn v 



V 2 dV = 4irp(Y)V 2 dV 



(4) 



V by definition is the velocity field smoothed by the win- 
dow function, once the smoothing scale is sufficiently large 
the distribution of V shall be very close to Gaussian so 
that Maxwellian dis tribution could be invoked to model p(V) 
(IBahcall etal.l[l99l 

MVW=/f(^) 3/ V 2 exp(-g)dy, (5) 



where the variance of V can be obtained with P vv , the power 
spectrum of v, through 



1 



(2tt) 3 



P vv W 2 d 3 /fc 



(6) 



Given that radius of the window function W(R) deployed to 
measure bulk flow is fairly large, one would expect that the 
velocity field smoothed at such scale will be well described 
by linear evolution of the initial condition, then if the initial 
distribution of velocity is Gaussian, for instance in the case 
of our simulation, Eq. [5] shall be a good approximation to the 
PDF of amplitudes of bulk flows. With the model, the most 
likely amplitude of bulk flow is simply V p = y / 2~/3oy, which 
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ranges of variance corresponding to different levels can be 
computed by the integral 

(/i)-/i«p(-g). m 

If v is curl free, a velocity potential field can be defined 
as 9{r) = -V -\/(Ha8) with the scale factor a = (1 and 

the growth rate B = d log D/d log a « QjJ/ 7 + ( 1 + m /2) Vt A /70 
(Hamilton 2001), we have 

The above expression relies on the assumption of negligible 
rotational velocity. In the linear regime, if the biasing of halo 
velocity to dark matter velocity is unity, it can be further sim- 
plified with the approximation Pee « Pi) where Pig is the 
linear matter power spectrum. 

3. BULK FLOW SHOWN IN THE PANGU SIMULATION 

3.1. The Pangu simulation and its halo catalogue 

The Pangu simulation (PS-I) is a large volume and high res- 
olution simulation, carried out under the scheme of the Com- 
putational Cosmology Consortium of China (dubbed C4). PS- 
I assumes a A cold dark matter (ACDM) cosmology model 
with parameters 

n m =0.26, fii = 0.044, n A = 0.74, 
/i = 0.71, cr 8 = 0.8, n s =l . 

The simulation contains dark matter only, and uses N=3, 072 3 
particles to follow the distribution and evolution of dark mat- 
ter within a periodic box with L = 1000/r'Mpc on a side. Each 
particle has a mass of 2.48915 x 10 9 /! _1 M©. The Plummer- 
equivalent force softening length is kept constant of 7/i _1 kpc. 

The PS-I starts from redshift Zmit = 127, initial positions 
and velocities of particles is generated with Zel'dovich ap- 
proximation from a glass-like particle set. Input linear power 
spectrum is computed with the CAMB dLewis et al.l2000l) . The 
simulation is then run with L-GADGET , a memory-optimized 
version of GAD GET 2 dSpringell 120051) . L-GADGET2 is de- 
signed to meet requirements of high performance computa- 
tions, only the tree-particle mesh algorithm is included to cal- 
culate the gravitational forces efficiently. Totally 64 snap- 
shots are saved from zm = 127 to redshift z = 0. The PS-I is 
performed on the supercomputer Lenovo Deepcomp7 
at Supercomputing Center of Chinese Academy of Sciences. 
We use 2048 cores and about 6TB memory at the peak time. 
The simulation consumes approximately 6.5 x 10 5 CPU hours 
(about 13 days) in total and consists of 6151 time steps. 

Dark matter halos are identified on the fly during the simu- 
lation for each snapshot, using the standard friends-of-friends 
(FOF) algorithm with linking length of 0.2 times the mean 
particle separation. Each FOF group must contain at least 20 
particles, at redshift z = there are 4.2 x 10 7 identified par- 
ticle groups. There are a little bit of ambiguity in definition 
of halo mass. Halo mass mostly used in literature is defined 
as the mass enclosed by a sphere centered on the halo center 
with certain radius, within which the average density is some 
factor larger than the critical density. A more convenient def- 
inition of halo mass is the total mass of all member particles 
of the FoF group, which is used in this paper. For simplicity 



the average of velocities of member particles are taken as the 
velocity of the halo. 

In this work we used a subset of the group catalogue as 
our full halo catalogue, which consists of FoF groups of 
mass larger than 2.5 x l0 n h~ 1 M@ only and makes 1.28 x 10 7 
entries in the end. Selected FoF groups contains at least 
100 member particles, corresponding Poisson fluctuation is 
greater than 10 which is the threshold used for FoF group 
identification, limiting discreteness error to 10% level. An- 
other reason is that it is not easy to detect large number of 
low mass halos in observation. Our prudence test verifies that 
including FoF groups with number of particle < 100 does not 
make any significant amendment to final results, even though 
effects of non-linearity are stronger. More importantly we will 
show that statistics of bulk flow a re no t sensitive to the mass 
of halos in the sample (subsection l3.41 i. 

3.2. The probability distribution function of bulk flow 

To ensure fair sa mpling to the simulation, the over- 
sampling algorithm of Szapudi) (11998b is implemented to gen- 
erate ~ 10 6 cells for R within 25- 100/r'Mpc, bulk flows of 
halos in random cells are estimated with Eq. Q] In our mea- 
surements we adopt mainly the spherical top-hat window, i.e. 
w = 1 for all halos inside window and w = otherwise, mean- 
while shell window function is also deployed for consistency 
check. Adopting the type of top-hat window function is just 
to simplify computation, in principle one could try Gaussian 
or other more sophisticated window functions, but it will not 
introduce change to the main results. Since our simulation 
box is limited in a cubic box of side length l/r'Gpc, probe of 
bulk flows in cells of radius > 100/i~'Mpc would be statisti- 
cally unreliable. 

Our measurements of bulk flows are displayed in Figure Q] 
For top-hat window, ay decreases with cell radius, the possi- 
bility of finding extremely large speed of bulk flow becomes 
smaller for larger volume. For shell window defined by two 
radius Rq < Ri, if Ro is not very close to R\, av(Ro,Ri) ~ 
oy(J?i). 

More importantly, PDFs of amplitude of bulk flow in sim- 
ulation, no matter measured with spherical top-hat window 
or spherical shell window, are all well described by the linear 
model. Agreement between simulation and model is better for 
larger volume as expected, linear theory slightly over-predicts 
try since nonlin ear Pee is lower than Pig at k ~ 0.1/iMpc _I by 
~ 20% already (Scoccimarro 2004). The comparison clearly 
lead to the conclusion that to a good precision p(V) obeys 
Maxwellian distribution which is completely determined by 
ay, what really matters is not the exact shape of the widow 
function but rather the corresponding ay. This lays out the 
solid ground for us to put different kinds of measurement to- 
gether for comparison. 

3.3. Bulk flow as mass weighted average of halo velocities 
It is known that attenuation to CMB temperature resulted 
from kSZ effect AT^sz °c v • I J n e dl in which l is the unit vec- 
tor of the line-of-sight and n e is the density of free electrons 
in the galaxy cluster. If the aperture used to measure kSZ ef- 
fect is sufficiently large, and the number of hot electrons in 
the cluster can be taken for granted proportional to the mass 
of host halo m/,, the total kSZ effect induced temperature fluc- 
tuation will be proportional to m/,v • l, thus the bulk flow esti- 
mated via kSZ effect of galaxy clusters in fact is in principle 



4 




Figure 1. Left panel: p(V) measured in spherical top-hat windows with different radius, solid lines are measurements from N-body simulation, from left to right 
with decreasing heights corresponds to radius of window function R = 100,50, 25/i~'Mpc respectively, dashed lines are prediction of the model of Eq.|5]& [8] 
with Pgg = Pi). Middle panel: p(V) measured in one spherical shell window defined by two radius (Rq,Ri) = (25, 100)/r'Mpc. Right panel: comparison of ay 
predicted by the linear theory (solid line is of spherical top-hat window, dashed line is of spherical shell window with Rq = R, R[ = 100/i~'Mpc) with estimation 
from simulation (blue crosses), and the ay which provide the best fitting to PDFs of bulk flow in simulation with Eq.[5](red squares). 
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Figure 2. Differences between V and V m in simulation. R is the radius of the spherical top-hat window function. Left: the relative difference between p(V,„) and 
p(V), the large fluctuation at the large V is due to the almost zero values of PDFs at tail. Middle: PDFs of amplitude difference between V and Vm in the same 
cell. Right: distribution functions of the angle between V and Y m in the same cell, a = cos -1 [V • Vm/(Wm)]. 

volume. If the sample volume is small, it could appear that V,„ 
deviates from V by ~ lOOkms" 1 though the possibility is tiny. 
We also notice that V m - V does not show apparent trend with 
V or V m . Pointing of mass weighted bulk flow does not coin- 
cide with V. The most likely angle between them is around 
1-3 degrees for top-hat window of R G (25, 100)/i _1 Mpc, and 
becomes smaller for larger volume. Note that the distribution 
of the difference angle has a rather long tail, for instance if 
R = 50/r'Mpc the probability of misalignment greater than 
10° is - 7.3%, yet not trivial. 

Hitherto only the kSZ measurements can provide estima- 
tion of mass-weighted bulk flow, meanwhile mass weighting 
does not introduce significant statistical diff erences, which i s 
actually supported by real observation (lLavaux et al.ll2.Q12h . 
thus hereafter we will just concentrate on the unweighted bulk 
flow. 

3.4. Halo mass dependence 

There is the possibility that bulk flow may depends on the 
typical mass of halo sample. The full halo catalogue is then 
divided into six subsamples by halo mass, measured oy is 
plotted in Figure [3] as function of the mean halo mass of the 
subsample. If the smoothing scale is large enough there is lit- 
tle dependence on mass of sampled halo of bulk flow, but if 
the window is too small to suppress non-linearity, e.g. the case 
of R = 25/r'Mpc, subsamples of low mass halos demonstrate 
slightly lower o v than high mass halo subsamples, indicating 
that low mass halos suffer of more nonlinear action than high 



V m is ratio of two Gaussian random variables, the total mass 
M = y^.OT,- and momentum P = ^ m{Vi, with the results of 
iPham-Gia et al.l (l2006h it is possible to work out a linear the- 
oretical model for p(V m ). Exact calculation needs knowledge 
of power spectra of matter, momentum and the correlation 
function between matter and momentum, which is not an easy 
task. However a quick inspection could give us a rough pro- 
file. V m can be expressed as (P / (M))/(M / (M)), i.e. numer- 
ator and denominator both are normalized by the mean mass 
enclosed by the window function, the variance of P/ (M) is 

dominated by the a v = ( 27r)~ 3 f P^W^d^k i f the smoothing 
scale is sufficiently large (Park & Park 2006), while variance 

of M/(M) is determined by a] = (2tt)" 3 / P ss W 2 d 3 k. Since 
(Jy 3> <Jg, p(M/(M)) behaves like a Dirac <5-function com- 
pared with p(P/ (M)), M can then be deemed roughly as a 
constant and p(V m ) f=s (M)p(P). P is a Gaussian random vari- 
able with variance a\ ~ (M) 2 (J V in linear regime, therefore 
p(P) ~ p(X)/(M) and then p(\ m ) - p(V). Our results of 
simulation data indeed reveal that there is little difference be- 
tween p(V m ) and p(V) (Fig. [2j. 

But in an individual cell V m does differ from V, both in di- 
rection and amplitude. As we can see in Fig. [2] p(V m — V) has 
width of several tens kms" 1 which decreases with larger cell 
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Figure 3. Halo mass dependence of ay . Symbols are measurement, their x- 
axis coordinates are the mean halo mass of subsamples. Lines are prediction 
of linear theory, below which are labelled with corresponding cell radius. 



mass halos. 

However for an individual volume, diversity in bulk flows 
given in different halo mass bins is very likely to be apparent 
when V <ay, although it declines with increasing V and turns 
out insignificant for V > 2oy. Our experiments disclose that 
differences between bulk flows of the lightest halo subsam- 
ple (thf, = 3.015 x lQ n h~ l Mo) and the mostly massive halo 
subsample (thf, = 1.374 x 10 i3 /i _i Mq) in numerical figures is 
actually extremely similar to Figure |2] Therefore if sampled 
volume in observation has low depth, for instance the effective 
R < 50/z _1 Mpc, and measured bulk flow is not of extremely 
high amplitude, attempt to find rigorous convergence among 
bulk flows estimated from different galaxy or galaxy cluster 
samples could be probably fruitless. 

3.5. Consistency between observation and model 

Systematical biases in dark flow, the bulk flow measured 
at high redshift, are not fully understood and precisely con- 
trolled, so we refrain ourselves from discussing high redshift 
case. Most of the local (or nearby) bulk flow measurements 
has redshift less than ~ 0.06, resulting redshift evolution of 
ay with respect to z = is of magnitude of a few percents at 
most, which can be comfortably ignored. Window functions 
in different works are not the same at all, but the excellent 
performance of the linear model provides an unified scheme. 
Since PDFs of bulk flow is solely determined by cry, indepen- 
dent of the type of the window, the radius of a top-hat window 
which gives the same linear oy as the window function used 
in observation can acts as the effective scale corresponding to 
a particular sample. 

To check if an observed bulk flow is consistent with ACDM 
model, we need to figure out the variance ranges of V. 
The most likely amplitude of V is V p = \f2~j3ay, derived 
via dp(V)/dV = 0, and the variance range at different lev- 
els are computed through Eq. [7] The Maxwellian distribu- 
tion is highly skewed, the one side probability /„ " p(V)dV = 
42.8%, but confidence probabilities at levels of la, 2a and 
3a are ei.2,3 = 68.3%, 95.5%, 99.7% respectively, so variance 
ranges at 2a and 3a significance are defined by (0, V2.3) with 

Jo 23 ' p(V)dV = £2,3, while the la variance is within (V{,Vi) 
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Figure 4. Measured local bulk flows against prediction of linear model. 
Symbols are some recently estimated amplitudes of local bulk flows from 
observational data of galaxies (Feldman et al. 2010; Nusser & Davis 2011) 
and supernoy ae I Haugbolle et al. 20071 IColin et aljhoi It IDai et alj|2ul*lt 
Wevant et al. 201 1). Solid line is the model predicted most likely bulk speed, 
dashed lines indicate variance at la level (68.3%), dot-dashed line and dotted 
line are at levels of 2a (95.5%) and 3a (99.7%) respectively. 

given J^p(V)dV = f£ p(Y)dV = e i /2. 

Several recent measurements from observational data are 
over-plotted upon model prediction in Fig. [4] One has to keep 
in mind that our calculation of effective scales has assumed 
that all other factors affecting estimation have been perfectly 
corrected, such as non-linearity, selection function and sky 
coverage incompleteness. Construction details of these sam- 
ples are often too sketchy to render appropriate weights for 
our calculation. However the imperfection actually reduces 
their effective volumes so that the true effective scales will be 
smaller, i.e. the data points will shift leftward along horizontal 
axis in Fig. [4] So we shall deem Fig.|4]as the mostly conserva- 
tive judgment of the consistency between observation and the- 
ory. Nonetheless, from Figure[4] it appears that observation of 
our local Universe does no t rule o u t the ACDM model. Those 
results of|H augb0lle "eTaT] (120071) . iFeldman et"aT] (120101) and 
IWeyant et al.1 (1201 ll) that often quoted as supporting evidence 
disfavoring standard ACDM model are around 3a level, but 
the significance will be smaller if error bars are taken into ac- 
count. In addition, considering that many other m easurements 
(including the not officially published report of Wand 120071) 
are in fact consistent with ACDM model, we prefer to choose 
conservative standpoint on the issue. 

3.6. Bulk flow and mass distribution in the cell 

It is interesting to investigate the relation between bulk flow 
and the mass distribution in the sample volume, one might 
wonder whether one could infer bulk flow from mass distribu- 
tion if peculiar velocity data is absent, since in linear regime 
Fourier modes of velocity field can be derived from modes 
of the density field. However from Eq.|2]it is clear that bulk 
flow is determined by those modes of wavelengths larger than 
the characteristic scale of the window function, if the Fourier 
transformation of the density field is restricted to the same 
volume in which bulk flow is measured, those modes of long 
wavelength accounted for bulk flow are missing. 

Our measurements confirm the expectation. The first quan- 
tity we checked is the mass monopole, the total mass (J2i m d 
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Table 1 

Correlation coefficients between amplitude of bulk flow and mass 
distribution 



Cell radius 


Mass monopole 


Mass dipole 


R(/r'Mpc ) 


£,m, N 


E'"; r i/'7 /r,\ 

Em W;/ 


25 


0.024 -0.0004 . 


0.139 0.134 


50 


0.004 -0.0188 . 


0.188 0.187 


100 


0.014 -0.0031 . 


0.245 0.249 



0.014 




angle between bulk flow and mass dipole (deg.) 

Figure 5. Probability distribution function of the angle between bulk flow 
and mass dipole. Dash-dot line is the expectation of null correlation. 

or the total number (AO of halos in the volume, which is equiv- 
alent to the density fluctuation smoothed by the window func- 
tion. Correlation coefficients are computed to denote the cor- 
relation strength between amplitudes of bulk flow and mass 
monopole (Table [TJ. Apparently bulk flow is not correlated 
with the total mass and the total number of halos in the vol- 
ume at all. 

Two kinds of dipoles of halo distribution are measured, 
(r,/r,) and the mass weighted one (^m,r,/V,)/^m,-. Our 
results show that including halo mass or not makes little dif- 
ference. As can be seen in Table Q] and Figure [5] mass dipole 
correlates with bulk flow very weakly both in amplitude and 
direction. The correlation becomes slightly tighter as sam- 
ple volume increases, the peak of the distribution function of 
the misalignment angle between mass dipole and bulk veloc- 
ity is at ~ 45° for R = 25/T 1 Mpc while shifts to ~ 32° for 
R= lOO/r'Mpc (FigureEJ. 

4. THE FASTEST BULK FLOW 

There are so me works claim dete c tion of unusually la rge 
bulk flow (e.g. iFeldmanet aDl20lot IWeyant et al.H201lT) . it 
is interesting to check properties of these special parts in 
ACDM universe. The largest amplitudes of bulk motion 
measured in our simulation for R = 25,50, 100/i -I Mpc are 
1070,778,514kms _1 respectively, the one for R = lOO/r'Mpc 
is already very close to those observational results. The pos- 
sibility of residing in the cell having the fastest bulk motion 
is defined by the ratio of the cell volume to the total vol- 
ume of simulation, m 0.42% for R = 100/r'Mpc. However, 
given the diameter of the cell as large as 200/r'Mpc against 



2.0 




0.5 - 



0.0 ' 

12 13 14 15 

L °gio m h (M Q /h) 

Figure 6. Difference in halo abundance between the fastest cell with R = 
100/j~'Mpc and the full halo catalogue, dridti/driAll, ^ n Cell * s the mass func- 
tion of halos in the fastest cell, dn A u is of the full halo catalogue. Red solid 
circles connected by solid line is of the fast cell, which is enclosed by dashed 
lines marking its Poisson variance. The blue dotted line is the average of the 
top-ten fastest cells. The horizontal solid line is of Anc e u / An A u = 1. 

the simulation box length l/r'Gpc, a huge volume moving at 
speed more than ~ 500kms _1 will yield observable features 
too prominent to be missed. 

In this report we will choose the R = 100/r'Mpc case as 
example to study peculiarity of the cell showing fastest bulk 
motion. The first physical quantity checked is the halo mass 
function in the cell, for comparison halo mass functions of the 
top-ten fastest cells (centers separated by at least 100/i _1 Mpc) 
are also measured. Shown in Figure [6] are measured cell halo 
mass functions divided by the halo abundance in the full cat- 
alogue. Among the ten cell mass functions, most (more than 
7 of 10) are smaller than the halo mass function of the full 
catalogue for m/, <^ IO^/t'Mq) In high mass regime due to 
the very small number of high mass halos we can not with- 
draw any reliable conclusions though the mass function of the 
fastest cells demonstrates a high tail. So far we would only 
cautiously conclude that in high bulk flow regions there is the 
tendency of finding less number of small mass halos. 

Distribution of all halo velocities in the cell is illustrated 
in Figure |7l also shown are velocity distribution functions of 
these halos in mass bins (0.25-0.3) x IO'V'Mq, (0.62- 
1.12) x \Q n h- l M Q and > 1.62 x 1O 12 /i _1 M . Mass binned 
halos do not exhibit any significant differences in aspect of 
velocity distribution, which eases the worry of possible bias 
in mass selected halo samples. Distribution function of the 
angle between halo velocity and bulk flow is very skewed to- 
ward small misalignment, the peak is around 20 ~ 30° but 
not the 0°, about 90% halos are moving in direction within 
60° to the bulk flow. It appears that halos in the cell with 
largest bulk velocity are more likely to have higher speed, the 
peak of velocity amplitude distribution of halos in the cell is 
at around 600km s" 1 while that of the full halo catalogue is at 
~ SSOkms" 1 (FigureEl). 

It has been examined that bulk flow basically is weakly cor- 
related with the internal mass dipole on average. But for the 
cell with the largest bulk flow, intuitively one would conjec- 
ture there should be certain very massive clumps neighboring 
to the cell, their gravitational action may play a dominant roll 
in causing such extreme bulk flow of nearby halos. As an at- 
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0.035 




Figure 7. Velocity distribution of halos in the cell with the largest bulk flow (R = lOO/i Mpc). Left panel is the distribution of amplitudes of halo velocities, 
think solid line is of halos in the cell, think solid line is of the full halo catalogue, other color lines are of halos in three different mass bins as labelled in figure 
legend. Right panel presents distributions of the angle between bulk flow V and halo velocity v/,, a = cos -1 [V ■ v/,/(Vv/,)]. 

Bulk flow and mass dipole I — • — • — — n • • • — ■ — ■ — • — n • ■ • — • — ■ — ■ 




Figure 8. Mollweide projection of directions of the largest bulk flow and 
mass dipoles. Mass dipoles are computed from the mass field represented 
by halos, ^m,r,/r,. The plus symbol marks the direction of the largest 
bulk flow (514kms-', cell radius R = 100/r'Mpc), the cross symbol is the 
direction of dipole of halos in the cell, the angle between them is 37.4°. Tri- 
angle, square, diamond, and circle symbols are directions of mass dipoles 
within shells (100, 150), (150,200), (200, 250), (250, 300)/r' Mpc, which de- 
viate from the bulk flow by about 30°, 33°, 21°, 17° respectively. The color 
map is the surface mass density contrast of halos in the cell, plotted in the 
Healpix scheme (dark color for high density contrast). 

tempt to justify the paradigm, mass dipoles in shells within 
R = 100-300/r'Mpc to center of the 100/i _1 Mpc cell with 
largest bulk flow are cal culated in four laye rs with the help 
of the Healpix package ( Gor ski et al.l 12005 ). if there is un- 
usual distribution of matter in a layer, mass dipole of the layer 
will be the efficient indicator. Projected directions of bulk 
flow and mass dipoles are displayed in Figure [8] When shell 
moves outward mass dipole pointing walks fairly randomly 
around the bulk flow, the misalignment angle varies between 
~ 20-40° which is analogous to the typical value in Fig- 
ure [5] It seems that dipoles of local environmental mass are 
only aligned crudely with the bulk flow. The correlation is 
not negligible, however since these misalignment angles are 
not small, we have no strong support from the simulation to 
attribute the extremely large bulk flow in such a huge volume 
mainly to inhomogeneous environment. 
Anomaly in clustering of halos in the special cell is detected 




o -0.1 I 

CJ 



1 10 too 

r (ir'Mpc) 

Figure 9. Correlation between bulk velocity and clustering of halos in the 
cell, cell radius is R = 100/r'Mpc. Top panel: solid line is the two-point cor- 
relation function £/, of halos in the cell with the largest bulk velocity, dotted 
line in red is the average over 1000 £/,s measured in randomly selected cells 
while the yellow shadow marks the corresponding lcr variance, dashed line 
in blue is the result of all halos in our halo catalogue. Bottom: correlation 
coefficients between amplitude of bulk velocity and t;i,(r) as function of r, 
summarized from measurements of 1000 randomly selected cells. 

(Figure |9j, halo two-point correlation functions of 1000 ran- 
domly located cells of radius lOO/r'Mpc and the full halo 
catalogue are also calculated. £/, averaged over the 1000 mea- 
surements agrees with £/, of the full halo catalogue at scales 
r < 30/r'Mpc, then drops down more quickly to zero at larger 
scales due t o the integral constrain t resulted from finite vol- 
ume of cell dLandv & Szalay 19931) . As we are interested in 
£/, in a finite volume, we did not bother ourselves to apply 
relevant correction. £/, of the cell with largest bulk velocity 
is higher than the average of random cells by around 2a at 
scales r < 20/i _1 Mpc, at larger scales the excess of cluster- 
ing power rises to level of ~ 2-3o\ In order to assess the 
statistical significance of the event, correlation coefficients of 
bulk flow amplitude and £/,(r) are computed from the 1000 
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random cells (bottom panel in Figure[9]l. Little correlation is 
detected at scales r < 20/i _1 Mpc, then at larger scales the cor- 
relation becomes much stronger, indicating that in cells with 
high bulk velocity it is truly more possible to find power ex- 
cess in halo clustering at large scales, which app ears to be in 
line with the findings in Macaulav et al.l (1201 11) . We further 
check the halo two-point correlation functions in the top-ten 
fastest cell, and we find that 7 of the 10 demonstrate power 
excess at similar scales. Note that the power excess at large 
scales can not be ascribed to integral constraint, for the lead- 
ing term of integral cons traint in this regime i s negative and 
will bend £/, downward (Landy & Szalayl ll993h . 

5. SUMMARY AND DISCUSSION 

Through analysis of the Pangu simulation, it is confirmed 
that bulk flow of halos follows Maxwellian distribution which 
is completely determined by a single parameter,the bulk ve- 
locity dispersion. We find that the dispersion measured in 
simulation agrees with the prediction of linear perturbation 
theory of structure formation very well, non-linearity only be- 
comes important when the sampling volume is very small. 
We find that the in most cases mass weighted bulk flow has 
some minor statistically differences to the unweighted one, 
but the will not affect the overall statistics significantly. It 
is also revealed that statistically bulk flow has little system- 
atical dependence on the mass of halos used for estimation. 
Based on the results, we propose a unified scheme to com- 
pare results from observational samples with theories. In the 
proposal, the scale at which bulk flow in a particular space of 
the Universe is measured is chosen to be the effective scale R 
which is the radius of a spherical top-hat window function W,i, 
that yields the same bulk velocity dispersion ay as the practi- 
cal window function Wo for the observational sample does in 
linear theory. Wo is not only determined by the sample geom- 
etry but also contains weights emerged from selection func- 
tion, incompleteness and etc., being analogous to the window 
function used in estimation of power spectrum. 

Variance ranges of bulk flow are clarified as well on the 
basis of Maxwellian distribution in the work, we make a 
rough comparison of some recently measurements with the 
ACDM model adopted in Pangu simulation, we find that part 
results do deviate from the model by about 3a but the tension 
between observation and model is not so strong as original 
works claimed. Estimated effective scales for observation re- 
sults in Figure |4] are in fact the upper limits, the true effective 
scales could be even smaller since we have assumed that those 
samples are of full sky coverage and their selection functions 
have been corrected for during estimation. Furthermore, ob- 
served bulk velocity consists of residuals from thermal motion 
of galaxies in their host halo, which is not included in calcula- 
tion of the velocity dispersion so far. More accurate modeling 
could be developed by assuming velocities of galaxies relative 
to their halos obey certain simple distribution, but it requires 
explicit knowledge of occupation details of galaxies in host 
halo, which in itself is already a challenging problem. A bet- 
ter way would be to deduct the random motion com ponent 
in the estimation procedure, such as the treatment in IWangl 
(l2007h . 

Correlation between bulk flow and dipole of internal mass 
is very weak, but is stronger for larger volume. If one happens 
to be living in a volume of radius greater than 100/r'Mpc with 
large bulk velocity, their observed mass dipole in the volume 
will have considerable chance of being unusually strong. Typ- 
ical misalignment angle between bulk flow and mass dipole is 



mostly likely around ~ 30-50°. This might introduce non- 
negligible systematical bias to cosmological probes involving 
local mass di stribution, such as the l ate-time integrated Sachs- 
Wolfe effect dRees & Sciamall 19681) . 

In our simulation there do exist volume of scale extending 
to 200/r'Mpc in diameter moving with extreme large bulk 
velocity more than 500kms _1 . Most halos inside the volume 
are moving in alignment with the bulk flow within 60°, and 
the flow shows no dependence on halo mass. Such group 
motion of numerous halos will generate prominent kSZ sig- 
nals, it is a rare event, but given its high speed and huge scale 
(200/r'Mpc versus l/r'Gpc), probability of detection is ac- 
tually not very small, which of course also depends on the 
inclination between bulk flow and line-of-sight. Another pos- 
sible observation effect is that galaxies in the volume could be 
dimmed or brightened on average by the extreme bulk flow 
than galaxies in other places, which is th e starting point of 
the eff ort tried bv lNusser et al.l d2011l) and [Abate & F eldmanl 
(HHI. 

Dipoles of mass outside the largest bulk flow region as en- 
vironment are not tightly aligned with the bulk velocity, but 
deviate from it by around ~ 20-40°. Interestingly we iden- 
tified that halo clustering of the particular volume is strength- 
ened apparently at scales r > 20/i _l Mpc, simulation results 
point out that such enhancement is not completely acciden- 
tal, at large scales two-point correlation function of halos in 
a finite volume is indeed mildly correlated with bulk velocity. 
Bulk velocity is dominated by Fourier modes of velocity at 
scales larger than the characteristic scale of the sample while 
in linear theory v(k) oc /k<5(k)/fc 2 , unusually large bulk veloc- 
ity seems to imply that there should be extraordinary super 
large mode of density fluctuation topping up in the region. 
However bulk flow is hardly correlated with the total number 
or mass of enclosed halos inside the sample volume, and as 
we checked the total number or mass of halos in the cell with 
largest bulk velocity is less than the mean value but still within 
2a variance range. Moreover, in linear regime Fourier modes 
of density fluctuation are independent, £ of halos inside the 
volume is controlled actually by modes of scale less than the 
characteristic scale of the sample. 

Aside from the theoretical puzzle, one question is whether 
the power excess of halo clustering in a region at scales 
> 20/r'Mpc can be used as indicator of candidate space of 
extremely large bulk flow, the advantage of using two-point 
correlation function is that clustering does not rely on di- 
rection of line-of-sight, the complication resulted from red- 
shift distortion in principle can be overcame by the ratio of 
£ h (r ~ 20-50/i _1 Mpc) to £ h at small scales e.g. ~ 10/r'Mpc 
where £/, is not correlated with bulk flow. A more serious con- 
cern is that if we are unluckily (or luckily) in a special region 
as large as our current largest galaxy survey with extremely 
large bulk flow, the measured clustering strength at and be- 
yond scale of baryonic acoustic oscillation would be signifi- 
cantly leveled up, could it be the case of power excess at very 
large scales in the Ba ryon Oscillation S pectroscopic Survey 
(BOSS) elaborated by lRoss et ail (I2012h ? To answer all these 
queries one surely needs multiple realizations of simulation 
of volume much bigger than our Pangu simulation, for the 
moment in this paper we have to leave these questions open. 

ACKNOWLEDGMENT 

This work is partly supported by the NSFC through grants 
of Nos. 10873027, 10873035, 11073055, and 11133003. 



Bulk flow in ACDM simulation 



9 



YPJ, WPL, XHY and PJZ are members of the Innovation 
group funded by NSFC (No. 1 1 121062). JP and XK acknowl- 
edge the One-Hundred-Talent fellowships of CAS. We appre- 
ciate stimulating discussion with Jiasheng Huang, Cheng Li, 
Guoliang Li and Lifan Wang. 

The Pangu simulation was carried out in the Supercomput- 
ing center of CNIC, CAS, under the collaboration scheme 
of the Computational Cosmology Consortium of China (C4), 
participating institutions are NAOC, PMO, SHAO and CNIC. 

REFERENCES 



Abate, A., & Feldman, H. A. 2012, MNRAS, 419, 3482 
Afshordi, N., Geshnizjani, G., & Khoury, J. 2009, JCAP, 8, 30 
Amanullah, R., et al. 2010, ApJ, 716, 712 
Bahcall, N. A., Cen, R., & Gramann, M. 1994, ApJ, 430, L13 
Colin, J., Mohayaee, R., Sarkar, S., & Shafieloo, A. 2011, MNRAS, 414, 
264 

Courtois, H. M., Tully, R. B., Makarov, D. I., Mitronova, S., Koribalski, B., 

Karachentsev, I. D., & Fisher, J. R. 2011, MNRAS, 414, 2005 
Dai, D.-C, Kinney, W. H., & Stojkovic, D. 201 1, JCAP, 4, 15 
Feldman, H. A., Watkins, R., & Hudson, M. J. 2010, MNRAS, 407, 2328 
Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, F. K., 

Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 759 
Hamilton, A. J. S. 2001, MNRAS, 322, 419 

Haugb0lle, T., Hannestad, S., Thomsen, B., Fynbo, J., Sollerman, J., & Jha, 

S. 2007, ApJ, 661, 650 
Jarosik, N, et al. 2011, ApJS, 192, 14 
Jha, S., Riess, A. G., & Kirshner, R. P. 2007, ApJ, 659, 122 
Kashlinsky, A., Atrio-Barandela, F, & Ebeling, H. 2011, ApJ, 732, 1 
Kashlinsky, A., Atrio-Barandela, F, Ebeling, H., Edge, A., & Kocevski, D. 

2010, ApJ, 712, L81 



Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 
Lavaux, G., Afshordi, N., & Hudson, M. J. 2012, arXiv: 1207. 1721 
Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473 
Macaulay, E., Feldman, H., Ferreira, P. G., Hudson, M. J., & Watkins, R. 

2011, MNRAS, 414, 621 
Mak, D. S. Y., Pierpaoli, E., & Osborne, S. J. 2011, ApJ, 736, 116 
Mersini-Houghton, L., & Holman, R. 2009, JCAP, 2, 6 
Moscardini, L., Branchini, E., Brunozzi, P. T., Borgani, S., Plionis, M., & 

Coles, P. 1996, MNRAS, 282, 384 
Nusser, A., Branchini, E., & Davis, M. 201 1, ApJ, 735, 77 
Nusser, A., & Davis, M. 201 1, ApJ, 736, 93 

Osborne, S. J., Mak, D. S. Y., Church, S. E., & Pierpaoli, E. 2011, ApJ, 737, 
98 

Park, C.-G., & Park, C. 2006, ApJ, 637, 1 

Pham-Gia, T., Turkkan, N, & Marchand, E. 2006, Communications in 

Statistics - Theory and Methods, 35, 1569 
Rees, M. J., & Sciama, D. W. 1968, Nature, 217, 511 
Ross, A. J., et al. 2012, arXiv: 1203.6499 

Sandage, A., Reindl, B., & Tammann, G. A. 2010, ApJ, 714, 1441 
Scoccimarro, R. 2004, Phys. Rev. D, 70, 083007 
Song, Y.-S., Sabiu, C. G., Kayo, I., & Nichol, R. C. 2011, JCAP, 5, 20 
Springel, V. 2005, MNRAS, 364, 1105 

Springob, C. M., Masters, K. L., Haynes, M. P., Giovanelli, R., & Marinoni, 

C. 2007, ApJS, 172, 599 
Strauss, M. A., & Willick, J. A. 1995, Phys. Rep., 261, 271 
Sunyaev, R. A., & Zeldovich, I. B. 1980, MNRAS, 190, 413 
Szapudi, I. 1998, ApJ, 497, 16 
Wang, L. 2007, arXiv:0705.0363 

Watkins, R., Feldman, H. A., & Hudson, M. J. 2009, MNRAS, 392, 743 
Weyant, A., Wood-Vasey, M., Wasserman, L., & Freeman, P. 201 1, ApJ, 
732, 65 

Wyman, M., & Khoury, J. 2010, Phys. Rev. D, 82, 044032 



